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ABSTRACT 

The inspirals and mergers of compact binaries are among the most promising events for ground- 
based gravitational wave (GW) observatories. The detection of electromagnetic (EM) signals from 
these sources would provide complementary information to the GW signal. It is therefore important 
to determine the ability of GW detectors to localize compact binaries on the sky, so that they can 
be matched to their EM counterparts. We use Markov Chain Monte Carlo techniques to study sky 
localization using networks of ground-based interferometers. Using a coherent-network analysis, we 
find that the Laser Interferometer Gravitational Wave Observatory (LIGO)-Virgo network can localize 
50% of their detected neutron star binaries to better than 50 deg^ with 95% confidence region. The 
addition of the Large Scale Cryogenic Gravitational Wave Telescope (LCGT) and LIGO- Australia 
improves this to 12deg . Using a more conservative coincident detection threshold, we find that 50% 
of detected neutron star binaries are localized to 13 deg^ using the LIGO-Virgo network, and to 3 deg^ 
using the LIGO-Virgo-LCGT-LIGO-Australia network. Our findings suggest that the coordination 
of GW observatories and EM facilities offers great promise. 



1. INTRODUCTION 

The era of gravitational-wave (GW) astron- 
omy is fast approaching. The advanced ver- 
sions of the Laser Interfer ometer Gravit a tional 
Wave Observatory (LIGO: : iBarish fc Weis^ ITOOl 
iSigg fc the LIGO Scientific Collaboration! 120081 ) and 
Virgo (jAccadiaet al\ |201H) are expected to make 
their first detections within the coming decade. 
Furthermore, construction has begun on the Large 
Scale Cr yogenic Gravitational Wave T el escop e 
(LCGT; iKuroda fc the LCGT Collaboration! [20101 ). 
and an additional advanced detector in Western Aus- 
tralia (referred to as LIGO-Austra ha or LAu in this 
work) is under serio us consideration (jBarriga et aZ.!!20l"ol : 
!Munch et aZ.![20Tl|) . Inspiraling and merging compact- 
object binaries, composed of neutron stars (NS) and/or 
stellar- mass black holes (BH), are promising sources 
for these detectors. For an advanced LIGO-Virgo 
network, predicted event rates for NS-NS binaries range 
from 0.4 to 400 ye ar"-'^ (with 40 bein g the "realistic" 
number given in (|Abadie et al\ \201C(] ) detectable to 
distances of several hundred Mpc, and approximately 
similar numbers apply for NS-10Mq BH binaries with 
detectable distances > I Gpc. Nearly face-on merging 
NS binaries are considered likely progenitors of shor t 
hard 7-ray bursts (SHBs; e.g., !EicIiler et al\ !1989! ). 
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!Metzger et"all ()2010[ ) suggest that NS bmary mergers, 
irrespective of their orientation, may produce radioac- 
tive decay powered transients with absolute magnitude 
peak luminosities of My ~ —15 in the optical at ^ 
day timescales. Thus, joint GW and electromagnetic 
(EM) observations can constrain the physics of generic 
NS binary mergers, the central engine and outflows 
of short GRBs, as well as cosmological parameters 
(see !Dalal et al\ !2006l: !Nissanke etall I2OIOI henceforth 
NIO). GW events should, therefore, be localized on 
the sky to sufficient accuracy to enable a match with 
their EM counterparts. Henceforth, we use the terms 
"sky localization" to refer to the measured sky position 
and its associated error uncertainty, while "sky errors" 
applies only to area estimates. 

For the typical NS binary inspiral event, whose du- 
ration is a few to tens of minutes in the detector fre- 
quency band, a single interferometer has a broad an- 
tenna pattern, and hence poor directional sensitivity. 
Two detectors restrict the sky localization to a single 
ring. In a network of three or more interferometers, rel- 
ative arrival times of signals at each detector allow for 
the reconstruction of the source's sky location to within 
a generally elliptical error area of 1-100 dcg^. For ini- 
tial a2d_advancedGW_d£^ networks, pre vious stud- 
ies (!Svlvestrel!2004l:!Cavalier e t al. 2006; Blair a/!!2008!: 
!Fairhurst!!2009l: !Klimenko et a/.n2011i ,Schutz.i201l!) have 
explored sky errors for sources with both unmodeled 
( "burst" ) and modeled waveforms using analytically de- 
rived timing formulae and Fisher-matrix methods. At 
high signal-to-noise ratio (SNR > 100) such methods 
are effective, providing lower bound error ellipses cen- 
tered on the true values of the parameters of interest. 
However, the majority of sources for LIGO will be de- 
tect ed near threshol d (SNR '--^ 8). As discussed in NIO 
and !Vallisnenl (pOOl . at "low" SNR (- 20-100), param- 
eter degeneracies can lead to quantitative and qualitative 
errors in the Fishcr-matrix approximation of the poste- 
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rior probability distribution functions (PDFs) for source 
parameters. Since early GW detections will likely be 
low-SNR, a full treatment of the sky localization will be 
critical when matching the GW event to an EM counter- 
part both (1) to ensure high-probability regions are not 
falsely excluded from consideration and (2) to minimize 
the number of false-positive transients that must be dealt 
with. Later, should event rates reach tens per year as pre- 
dicted, the ease of finding counterparts will be directly 
related to how well GW networks can localize events. In 
this paper we use Markov Chain Monte Carlo (MCMC) 
techniques to map the full PDF of the sky location where 
an event may have occurred. Specifically, we examine sky 
localization for non-spinning NS-NS populations for ad- 
vanced GW detector networks, including LIGO (the two 
4 km sites situated at Hanford and Livingston), Virgo, 
LCGT, and LIGO- Australia. For conciseness, we re- 
strict our analysis here to NS-NS binaries, results for 
NS-BH binaries being similar. For networks compris- 
ing the initial versi ons of LIGO and Virgo detectors , 
several works (e.g.. iRoyer. Mever fc GhristensenI 120071 : 
Ivan der Sluvs et aLll2009l : iRavmond et aZ.ll2009[ )lmve ad- 
dressed sky localization using MCMC techniques for sin- 
gle NS-NS and spin-precessing NS-BH systems. After 
the present article was submitted , a sim ilar MCMC study 
was performed by lAvlott et al\ (I2011D . They focus on 
the addition of LIGO Australia, and find consistent re- 
sults to ours. 

We follow similar techniques to those developed in NIC, 
which focuses on distance determination using advanced 
GW detector networks for sources with known EM coun- 
terparts. In contrast to NIO, here we instead ask how well 
a network of ground-based GW detectors can determine 
the sky positions of previously unknown compact binary 
sources, a central question when planning complemen- 
tary EM observations for counterparts. 

This paper is organized as follows. Section [5] out- 
lines the extraction of sky position from the GWs emit- 
ted by inspiraling binaries. Section |3] discusses our sky- 
localization results for individual systems and popula- 
tions of NS-NS binaries. Section |4] presents our conclu- 
sions. 

2. SKY LOCALIZATION DETERMINATION 

Ba, s ed on optimal matched filtering (lOppenheim et aTl 
[1981 lFmI][T99l ICutler fc FlanaganI iBglT " we extract 
the sky position n for each NS-NS binary using knowl- 
edge of the expected GW waveform, where n = (0, 0) 
is the vector pointing to a binary on the sky (the waves 
therefore propagate to the Earth along — n). The colat- 
itude 6 and longitude are related to the declination S 
and right ascension a, by0 = 7r/2 — (5 and 4> = a— CAST 
respectively, where CAST is Greenwich Apparent Side- 
real Time. We use only the early inspiral portion of the 
waveform, which for low mass syste ms provides most of 
the s ignal for advanced detectors ([Flanagan fc Hughei 
119980 and is modeled accurately using post-Newtonian 
(FN) expansions in general relativity. Specifically, we 
use the non-spinning restricted 2PN waveform in the fre- 
quency domain for the two GW polarizations and ; 
see Equations. (12)-(14) in NIO. The detector antenna 
functions depend on n and the binary's polarization an- 
gle. The overall amplitude of the GW waveform encodes 
the source's orientation, sky location, luminosity distance 



and redshifted chirp mass (see discussion in NIC). For 
simplicity we assume that the inspiral waveform ends 
abruptly prior to merger at the innermost stable circular 
orbit. The time of flight from source at direction n to 
detector at location r involves the scalar product n • r, 
and differences in time of flight among detectors in the 
network dominate sky localization. 

In order to infer the sky position (cos 0,0), we ex- 
plicitly map out the posterior PDF for all source 
parameters (including chirp mass, orientation, etc.) 
given an observed data stream at a detector, us- 
ing MCMC. The Metropolis-Hastings MCMC algo- 
rithm used is base d on a generic versi on of Cos- 
moMC, described in [Lewis fc Bridld ([2002[ ). For sim- 
plicity, we assume Gaussian, stationary, and zero-mean 
noise that is independent and uncorrclated between 
detectors. We take the anticipated noise sensitivity 
curve for a single advanced LIGO detect o r, giv en in 
IHarrv fc the LIGO Scientific Collaboration mm . for 
broadband tuning, to be representative of all our detec- 
tors, imposing a low-frequency cut off of 10 Hz. Since 
we trigger detections off the expected and not the ob- 
served SNR, we do not incorporate selection effects into 
our analysis. Therefore, we take prior distributions in 
all source parameters to be flat over the region of sam- 
ple space where the binary is detectable at an expected 
network SNR = 3.5. The expected network SNR is de- 
fined as the root sum square of the expected individual 
detector SNRs. For each MCMC simulation, we derive 
solid angle areas over (cos 0,0) for 68%, 95%, and 99% 
confidence regions. 

2.1. Detected Binary Populations 

The selection of our "detected" binaries is central to 
the derivation of represent ative sky loc alization statis- 
tics for binary populations (|Schutzll2011[ ) . We simulate a 
million binaries out to z = 1 , assuming a constant comov- 
ing v olume density in a ACDM universe ([Komatsu et al\ 
l2009t) . and with random binary sky positions and ori- 
entations. Each NS has a physical mass of IAMq. We 
consider two plausible scenarios where each binary in the 
total population is selected: (1) if its expected network 
SNR is greater than a network threshold of 8.5 (Case I) 
and (2) if each expected SNR at both LIGO Hanford and 
LIGO Livingston has a value of 6 or higher and if the ex- 
pected network SNR is above 12 (Case II). We choose 
to threshold off the LIGO detectors, as opposed to any 
other two detectors, only as a representative example out 
of all other possible coincident threshold scenarios. Case 
I uses a networ k thresholding criterion similar to that 
used in NIC and lSchutzl (|20Tlli . 

3. RESULTS AND DISCUSSION 

3.1. Individual Binaries 

We first examine sky localization for NS-NS binaries 
taken at random from our Case I or Case II detected 
samples. For each binary we assign a unique noise real- 
ization to each detector, which we keep constant when 
adding and subtracting detectors to a network. This en- 
ables a meaningful comparison between the performance 
of different networks. We consider networks comprising 
combinations of LIG O, Virgo, LAu, a nd L CGT. 

In a greement with lFairhurst[ ([20101 ) and iWen fc Cheiil 
(|2Q10|) . we find elliptically shaped errors for the majority 
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Fig. 1. — Sky localization for a high SNR binary under different 
network configurations, as labeled. The solid black curves indicate 
the 68% and 95% confidence regions (c.r.). Additional detectors 
increase the network SNR and decrease the error ellipse, however 
note the significant localization enhancement provided by LAu in 
particular. The origin (0,0) of each plot represents the source's true 
position, and the solid black lines denote the confidence regions. 

of our examined NS-NS binaries. This is unsurprising 
given that sky localization reconstruction is dominated 
by differing GW arrival times at each detector rather 
than the direction-dependent antenna functions in the 
GW amplitude. For the three-element LIGO-Virgo net- 
work, sources located toward the detector plane produce 
elongated error ellipses and have relatively poor angular 
resolution. However, the error ellipse significantly de- 
creases in size as additional detectors are added to the 
network. The inclusion of LAu is particularly favorable 
and breaks the LIGO-Virgo(~LCGT) degenerate plane. 
As a representative example. Figure [T] shows a relatively 
high SNR GW signal for a NS-NS binary located at 180 
Mpc with an inclination angle of cost = 0.7 and a sky 
location of (cos6l = -0.3, (/) = 2.9). The expected SNRs 
at LIGO Hanford; LIGO Livingston, Virgo, LAu, and 
LCGT are 6.7, 7.8, 12.4, 10.5, and 8.9, respectively. The 
specific orientation and shape of the ellipse are depen- 
dent upon the sky position and orientation, and noise 
realization. 

However, for a handful of events near threshold, we 
find that sky error areas are considerably non-ellipsoidal 
and exhibit multimodal distributions, in particular with 
networks that do not include LAu. As expected, their 
errors are not centered on the source's true sky posi- 
tion and we find that 5% of our selected binaries have 
true positions th at lie outsi d e the 9 5% confidence region . 
Analvses suc h as iFairhurstl (|2010( ). iWen fc ChenI (|2010f ). 
and ISchutd (|2011[ ) cannot reproduce such features due 
to limiting assumptions implicit in timing and Fisher 
information methods. Figure [2] shows an example of a 
low-SNR NS-NS binary located at 567 Mpc with an in- 
clination angle of cost = —0.93 and a sky location of 
(cos6i = -0.36,0 = 1.5). The expected SNRs at LIGO 
Hanford, LIGO Livingston, Virgo, LAu, and LCGT are 
5.4, 6.2, 3.1, 5.6, and 3.0, respectively. 

We find that error areas typically decrease by a factor 
of 4-5 by including LAu in the network. For particu- 
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Fig. 2. — Sky localization for a low SNR binary. Contour levels 
are as in Figurc[T] In this example, the three-detector network finds 
multiple, widely separated islands of high likelihood, offset from the 
true location. Additional detectors, however, can break underly- 
ing degeneracies. Furthermore, the binary's sky error at 95% c.r. 
increases slightly with the five-detector network compared to the 
four-detector (LAu) network, because in this particular instance, 
the noise realization is unfavorable at LCGT. 



larly weak signals, multimodal peaks and statistical bi- 
ases occur because multiple likelihood peaks for different 
values of the GW waveform's amplitude become indis- 
tinguishable from each other due to the uncertainties in 
the signal's time of arrival at each detector. 

3.2. Detected Samples of Binaries 

We now examine cumulative distribution in sky errors 
for ensembles of GW NS-NS events using different de- 
tector networks. We randomly choose events from our 
samples of detected NS-NS binaries usi ng th e two differ- 
ent selection criteria detailed in Section 12.11 

3.2.1. Case I: triggering on a Network of Detectors 

In Case I we set a total GW detector network thresh- 
old of 8.5, which implies an approximate SNR threshold 
per detector of 8.5 /Vb ^ 3.8 for each of the five detec- 
tors. As a first approximation, we estimate a detectable 
range of events within the maximum network capabil- 
ity: each detector in the five-detector network has a 
"weighted geometric average" (optimal) range of about 
420 Mpc (940 Mpc) for NS-NS eventsQ The geometric 
average statistic used here is a weighted angular average 
over all sky positions and orientations, which is a factor 
of '--^ 2.24 smaller than the optimal range for a face-on 
binary that is located directly above the detector (see 
iFinn fc ChmIol[l99l . 

In order to obta in detection event ra tes, we follow the 
approach given in lAbadie et al\ ()2010[ ) where detection 
ranges are derived by thresholding off a single LIGO 
SNR of 8. As the noise in r e ality is non-Gaussian and 
non-stationary, lAbadie et al\ ()201Cl[ ) use the range for a 



Similar to lAbadie eFai] | |20TOD . we do not henceforth incor- 
porate cosmological redshifts for our NS masses when estimating 
detectable ranges and rates. 
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single detector to represent the network of LIGO- Virgo 
detectors in order to achieve desired false alarm rates. 
For comparison purposes, we apply the same argument 
to our analysis: a five-detector network will have an ap- 
proximate threshold per detector of (8.5 x a/3) /\/5 ~ 6.6, 
and hence a geometric average (optimal) r ange of about 
240 M pc (540 Mpc) for NS-NS events. lAbadie et al\ 
(|2010[ ) provide a simple prescription to compute the GW 
NS-NS detection rate using NS-NS coalescence rates per 
galaxy, and the number of galaxies accessible within a 
given GW reach. The NS-NS coalescence rate per galaxy 
is estimated either by extrapolating from the observed 
sample of NS-NS binaries detected via pulsar measure- 
ments, or by using population-synthesis methods. Using 
Table II with the associated low-realistic-high~maximum 
NS-NS coa lescence rates per ga lax^fl and Equations (1) 
and (5) in lAbadie etai} (|2010f l. we define our realistic 
detection rate of ~ 65 NS~NS binaries per year. Cor- 
responding low, high and maximum detection rates are, 
respectively, 0.7, ^ 660 and ^ 2650 NS-NS binaries per 
year seen by the LIGO+Virgo+LAu+LCGT network. 

Turning to our results, we take a sample of 98 NS-NS 
binaries detected by the full five detector network, cor- 
responding to an observation time of ^ 8 months based 
on our realistic rate estimates. In Figure [3J we show 
the normalized cumulative distribution of sky errors in 
square degrees for our sample. We show distributions for 
subsets of systems detected by different networks which 
are normalized to the full sample illustrating the reduced 
number of detections. Table [1] illustrates sky errors for 
25%, 50% and 75% of NS-NS binaries from the sample 
detected by a particular network. Among notable fea- 
tures: (1) the addition of detectors to the network, in 
particular LAu, significantly reduces sky localization er- 
rors. We find that 50% of all detectable NS-NSs are 
localized at 95% confidence region to within 10-20 deg^ 
with any four or five detector network including LAu, 
and to within 110 deg^ with only the three LIGO-I- Virgo 
network. (2) The number of detected binaries doubles 
as the numbers of detectors in a network increases from 
three to five. 

3.2.2. Case II: triggering on Individual LIGO Detectors 

For the Case II scenario we define detection using the 
more stringent requirement of SNR > 6 at each LIGO 
detector. Using an SNR of 10.4 (~^ 6 x a/3) in each 
LIGO detector (whe re the factor of a/3 follows from the 
lAbadie et al\ (|2010[ ) correction for non-stationary and 
non-Gaussian noise), we compute a weighted geometric 
average (optimal) range of 150 Mpc (340 Mpc). Invoking 
a similar argument as in Case I, we estimate a realistic 
detection rate of ^ 17 NS-NS binaries per year. Cor- 
responding low, high and maximum detection rates are, 
respectively, 0.2, ~ 165, and ~ 660 NS-NS binaries per 
year. 

Similar to Case I, we take a detected sample of 88 
NS-NS binaries seen by the fuh LIGO-Virgo-LAu-LCGT 
network. Figure [4] shows the normalized cumulative dis- 

® lAbadie et al.\ 1)20101 ) assign rate estimates to one of four cate- 
gories as detailed in their Section IV; when rate PDFs are available, 
realistic refers to the mean of the PDF, low and high denote the 
95% pessimistic and optimistic confidence intervals, and maximum 
is the upper limit quoted in published literature. 
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Fig. 3. — Normalized cumulative distributions as a function of 
the sky-error area (square degrees) of a sample of NS-NS binaries in 
Case I detection scenario. Key: Solid/dotted lines denote 68% and 
95% c.r. respectively. Black: LIGO+Virgo+LAu-f-LCGT network, 
green: LIGO-|-Virgo-|-LAu, red: LlGO+Virgo-f-LCGT, and blue: 
LIGO+Virgo only. 

tribution of our NS-NS binary sample as a function of sky 
error in square degrees. We find that 50% of NS-NSs are 
detected at 95% confidence level to within 5 deg^ with 
any network including LAu, and ^ 15 deg^ with only 
the LIGO-Virgo network. As before. Table [T] shows the 
cumulative distribution of sky errors for NS-NS binaries 
detected by different networks. Despite differences in our 
analyses and astrophysical populatio n mod els us ed, our 
result s are consistent with iFairhurstI (|2010[ ) and ISchiit^ 
(I20T1I) . Once again, the inclusion of LAu substantially 
improves localization errors and the binary detection rate 
rises with more detectors. In contrast to Case I, however, 
we find that the number of detected binaries increases 
only by a factor of 1.3 going from a three- to five-detector 
network, because both LIGOs must have SNRs > 6. Us- 
ing all five detectors, the Case II scenario results in a 
factor of ~ 5 fewer binaries than in the Case I counter- 
part. We also find that the addition of a fourth or a fifth 
detector does not improve sky coverage as significantly 
in Case II as in Case I. 

4. CONCLUSIONS 

This work examines sky localization for both individual 
systems and populations of NS-NS binaries using differ- 
ent networks of advanced GW detectors. For the major- 
ity of optimally oriented NS-NS binaries examined, we 
show good agreement of our MCMC derived errors with 
the error ellipses obtained from analytical timing accu- 
racy and Fisher matrix formulae. However, for a handful 
of SNR signals at threshold (in particular with the stan- 
dard geometrically degenerate LIGO-Virgo network), we 
show that sky error regions can be non-ellipsoidal and 
non-contiguous, show multimodal distributions, and the 
best-fits can be shifted away from their true values. Of 
particular relevance, our results show that the inclusion 
of LIGO- Australia in a worldwide GW detector network 
improves localization errors both for individual and pop- 
ulations of binaries up to a factor of ^ 5, reducing the 
appearance of multimodal islands. Finally, the number 
of detected binaries increases with the number of detec- 
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TABLE 1 

Sky Errors in deg-^ at 68%/95%/99% Confidence Regions (c.r.) for NS-NS Binary Populations Detected 
Using Two Different Selection Criteria with Varying Networks. Percentages of detected binaries as a 

function of detector network are also indicated. 



Network 


Fraction of binaries 


LIGO+Virgo (LLV) 


LLV+LAu 


LLV+LCGT 


LLV+LAu+LCGT 






68%/95%/99% 
in deg^ 


68%/95%/99% 
in deg^ 


68%/95%/99% 
in deg^ 


68%/95%/99% 
in deg^ 




25% of binaries 


11/29/45 


2/5/10 


4/12/20 


2/7/10 


Case I 


Median percentile 


17/53/88 


4/13/22 


9/30/49 


4/12/20 




75% of binaries 


30/124/200 


7/30/53 


15/47/82 


7/22/39 




25% of binaries 


2/9/15 


1/2/3 


1/3/5 


1/2/3 


Case II 


Median percentile 


5/13/21 


2/5/9 


3/7/11 


1/3/5 




75% of binaries 


12/36/60 


3/9/14 


5/20/33 


2/7/12 


Case I 


% of detected binaries 


49 


66 


69 


100 


Case II 


% of detected binaries 


73 


97 


83 


100 



events result in error areas of ^ lOdeg^. In the optical, 
cross-correlating localization error areas > 10 deg^ with 
as complete as possible galaxy cata logs, such as the "Lo- 
cal U niverse" census proposed by iKulkarni &: Kasliwal 
(l2009f) . should aid EM search strategies ([Kasliwal et al. 
1201 Ih . As recent observations indicate, a small number 
five at present) of SHBs appear to be located several 
tens of kpc away from their host galaxies (e.g., iBergeii 
I2010t iFong et al. |[20Tof ). Thus, search strategies wiU 
need to include the possibility of such effects when look- 
ing for transients of NS binary mergers. Moreover, time- 
domain surveys provide an estimate of the false-positive 
rate of dynamic galactic (foreground) and extragalactic 
(background) transients. Currently operating and future 
optical and radio EM facilities are thus capable of prepar- 
ing for and performing follow-up of GW events. 
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